What is needed of a tachyon if it is to be the dark energy? 
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We study a dark energy scenario in the presence of a tachyon field <f> with potential V(<f>) and a 
barotropic perfect fluid. The cosmological dynamics crucially depends on the asymptotic behavior of 
the quantity A = — M P V^/V 3 ^ 2 . If A is a constant, which corresponds to an inverse square potential 
V{4>) oc </> -2 , there exists one stable critical point that gives an acceleration of the universe at late 
times. When A — > asymptotically, we can have a viable dark energy scenario in which the system 
approaches an "instantaneous" critical point that dynamically changes with A. If jA| approaches 
infinity asymptotically, the universe does not exhibit an acceleration at late times. In this case, 
however, we find an interesting possibility that a transient acceleration occurs in a regime where |A| 
is smaller than of order unity. 

PACS numbers: 98.70.Vc, 98.80.Cq 



I. INTRODUCTION 

There are few more taxing questions facing cosmology 
today than what is the nature of the dark energy in the 
Universe? This almost uniform distribution of energy 
density with a substantial negative pressure completely 
dominates all other forms of matter and yet the best we 
can do is to infer its existence from data (see Refs. pj for 
reviews). Moreover, because today it appears to behave 
just like a cosmological constant, we are naturally led to 
take seriously the possibility that what we are observing 
is a remnant of the physics of the very early Universe! 

Over the past few years there have been many papers 
devoted to addressing the nature of the dark energy. Ex- 
planations include: a true cosmological constant possibly 
arising from a Landscape type picture of string vacua [2j ; 
dynamical 'Quintessence' fields Q which have tracker like 
properties where the energy density in the fields track 
those of the background energy density before dominat- 
ing today; K-essence scenarios |4j , where the acceleration 
is driven by modified kinetic terms in the underlying ac- 
tion; modifications of gravity |{| (mainly motivated by 
Braneworld models) which lead to late time accelerat- 
ing solutions of the modified Friedmann equation and 
Chaplygin gas models |(| which attempt to incorporate a 
unified description of dark energy and dark matter. The 
list goes on, but the aim of all the models is the same, 
to explain why the Universe has only recently started ac- 
celerating with an energy density so close to the critical 
density. 

In this paper, we turn our attention to the issue of the 
tachyon as a source of the dark energy. The tachyon is 
an unstable field which has become important in string 
theory through its role in the Dirac-Born-Infeld (DBI) ac- 
tion which is used to describe the D-brane action 0, @ ■ 



A number of authors have already demonstrated that 
the tachyon could play a useful role in cosmology 0, 
independent of the fact that it is an unstable field. It 
can act as a source of dark matter and can lead to a 
period of inflation depending on the form of the asso- 
ciated potential. Indeed it has been proposed as the 
source of dark energy for a particular class of potentials 
ES IH EE However, there has not really been 
an effort to understand the general properties of tachy- 
onic cosmologies. We attempt to do that here. Starting 
with the four-dimensional Dirac-Born-Infeld (DBI) ac- 
tion, the tachyon field (with arbitrary potential V((f>)) 
is coupled to a background perfect fluid of radiation or 
matter. Modifyin g th e procedure introduced in [1 5[ (see 
also Refs. [H 111 111 E E E E E E EH! on 
a related theme), the evolution equations for the sys- 
tem of tachyon plus a background fluid can be written 
as first order differential equations involving two vari- 
ables x and y where x oc <j> and y oc y/V(<p) /H, H be- 
ing the Hubble parameter. Such equations have solu- 
tions whose behaviour depends crucially on the quantity 
A = -MpV^/V 3 / 2 , where = dV/d(j). We analyse this 
behaviour for a wide class of potentials. Depending on 
the asymptotic behaviour of A we either obtain late time 
attractor solutions which correspond to the well known 
inflationary cosmologies associated with the tachyon po- 
tential V{<f>) oc 1/0 2 , or perhaps more interesting cos- 
mologies which include (for A — > 0) asymptotic behaviour 
where the system approaches an "instantaneous" critical 
point that dynamically changes with A. For the case 
where |A| — > oo asymptotically, the system does not lead 
to late time acceleration. However, it does show a period 
of transient acceleration in the regime where |A| < 0(1). 
It raises the possibility that we are living in such a tran- 
sient regime. 
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The rest of the paper is as follows: Section II intro- 
duces the DBI action and the associated equations of 
motion for the tachyon field including the background 
perfect fluid. The equations of motion in terms of the 
new variables x, y and A are presented along with the 
effective equation of state for the tachyon. Potentials 
which demonstrate the particular asymptotic features of 
A = const, or oo are included. Section III considers 
the particular case of constant A finding all the associ- 
ated fixed points for the equations of motion and demon- 
strating their stability before going on to show how the 
existence restricts the allowed form of the potential for 
the tachyon. We show how easy it is to obtain late time 
accelerating solutions in this case. Section IV considers 
the case A — » for both a non-oscillating and oscillating 
late time evolution of 0, considering a class of potentials 
that have these features. As before we find it is possible 
to have late time acceleration and classify the conditions 
that have to be met by the tachyon potential. Section V 
does a similar thing but for a class of potentials which 
lead asymptotically to A — > oo. Although generally we 
do not have late time inflation, we show there is a novel 
feature present which is a transient period of inflation for 
exponential potentials. We summarise in Section VI. 



The background equations of motion are 



H = 



2M^/1 



1^- (3) 
2 Ml ' 1 ' 



l-(/) 2 

p B + 3-fHpB = 



3^+^ = 0, 



(4) 
(5) 



together with a constraint equation for the Hubble pa- 
rameter: 



3M 2 H 2 = 



1-4) 



PB 



(6) 



Let us rewrite the above equations in an autonomous 
form. We define the following dimensionless quantities: 



= Hcf>'. 



(7) 



where a prime denotes the derivative with respect to the 
number of e-folds, N = In a. Then we obtain the follow- 
ing equations: 



II. DBI MODEL 



The Dirac-Born-Infeld type effective 4-dimensional ac- 
tion for our system is described by 



x = -(1 - a; 2 )(3x- V3Ay), 
yi=y(-V3Xxy-^-^ 2 

2 V 

A' = -VSX 2 xy(T - 3/2) , 



3 7 



(8) 
(9) 
(10) 



where 



S = 



rf 4 :, 



M 2 
— v p 



-V(4>) \J- detig^ + d^d^ 



(1) 



where M p is the reduced Planck mass, R is the scalar 
curvature and V(<j>) is the potential of the tachyon field 
<j>. The above tachyon DBI action is believed to describe 
the physics of tachyon condensation for all values of <f> as 
long as string coupling and the second derivative of <f> are 
small. 

We shall consider a cosmological scenario in which the 
system is filled with the field tj> and a barotropic perfect 
fluid with an equation of state ps — {"f~l)pB- Note that 
7 = 1 for a pressureless dust and 7 = 4/3 for radiation. 
In a spatially flat Friedmann-Lemaitre- Robertson- Walker 
(FLRW) metric with a scale factor a, the pressure and 
the energy densities of the field <j) are given, respectively, 
by 



V{<t>) 



(2) 



A = - 



M p V^ 



V 3/2 

From Eqs. © and JfJJl we have 



r = 



w_ _ 

\\- also define 



7' 



(7 - x 2 )y 2 



n<4 = — — 

Per 



3M 2 H 2 \/1 



y/T^ 2 



Per 



PB 



3M 2 H 2 1 



(11) 

(12) 

(13) 
(14) 



where p CI — 3M 2 H 2 . Note that these satisfy the con- 
straint equation + Qb = 1 by Eq. ©. Since < 
Q$ < 1, the allowed range of x and y is < x 2 + y 4 < 1. 
Therefore both x and y are finite in the range < x 2 < 1 
and < y < 1. The effective equation of state for the 
field 6 is 



70 



12 



(15) 
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which means that 7^ > 0. The condition for inflation 
corresponds to (j) 2 < 2/3 It is also convenient to 
introduce a standard deceleration parameter as 

_ ' da 3 -, 3 7-x 2 2 
a 2 - 2 2 VI - x 2 

Equation (|10|) implies that one can have A = const for 
T = 3/2. In this case integrating Eq. (JH} gives (lflllljl 

= M 2 <j)- 2 . (17) 

This corresponds to the potential for scaling solutions in 
the context of braneworld cosmology [2lt l23l |27| . Re- 
cently a phase-space analysis was performed in Ref. 01 
for the same potential in the tachyon system. In this pa- 
per, we are intend to broaden this type of investigation to 
allow for the more general situation where the potential 
is not restricted to Eq. l|17fl. In fact the tachyon poten- 
tials we introduce are motivated either by string theory 
or based purely on phcnomcnological considerations. 

For example, in the case of a general inverse power-law 
potential [Tal^ 

V(4>) = M i ~ n <p~ n , (18) 

one has A cx 0("~ 2 )/ 2 and T = (n + l)/n. Therefore A is 
constant for n = 2, but it dynamically changes for k/2. 
In the limit of <fi — > 00 we have A — > for < n < 2 and 
A — > 00 for n > 2. 

It is convenient to classify the potential V(<f>) depend- 
ing on the asymptotic behavior of A as in the case of a 
normal scalar field 0. They can be classified as follows: 

• (i) A = const. 

In the tachyon system, the inverse square poten- 
tial V{4>) = M 2 <fi~ 2 gives a constant A. Note that 
an exponential potential corresponds to a constant 
value of A = — M P V^/V in the case of a normal 
scalar field [Til ]. 

• (ii) A — » asymptotically. 

There exist a number of potentials that exhibit this 
behavior. For example 

- (iia) V = M 4 - n <j)- n with < n< 2. 

In this case A (cx <ft( n ~ 2 ^ 2 ) approaches as 
(j) — > 00. It is known that inflation occurs for 
< n < 2 [13. 

- (iib) V = Voe 1 /^. 

This is the potential that is used in 
the quintessence scenario Q, and A cx 
e - 1/(2^0) j wmc ] 1 satisfies A — > as 4> — > 00 ■ 

- (iic) V(4>) = Vbel^V. 

This potential may appear as the excitation of 
massive scalar fields on the D-brane [2|j and it 
has a minimum at </> = 0. Since A cx (j>e~i M ^ 
in this case, A — > as — > 0. 



• (iii) |A| — > 00 asymptotically. 

There are also several potentials of interest that 
give this behavior 

- (iiia) V = M 4 - n (/)- n with n > 2. 

As seen earlier, in this case A (cx <ft( n ~ 2 ^ 2 ) 
approaches 00 as <f> — > 00. 

- (iiib) V = V e-^. 

This potential was considered in Ref. [3(j in 
the context of tachyon inflation. Unlike the 
case of conventional cosmologies, in tachyon 
cosmology, the exponential potential does not 
lead to a constant A. Since A is proportional 
to 65*"^, we find A — > 00 as — > 00. This po- 
tential may appear in the late time behaviour 
of D3 anti-D3 cosmology |3l|. The tachyon 
in the D3 anti-Da system is a complex scalar 
field T = (j)e ie fig. When 6 is constant, the 
effective action of coincident D3 anti-D3 is the 
same as Eq. |QJ 33j . The tachyon potential for 
the coincident D3 anti-D3 is the same as the 
tachyon potential of non-BPS D3-brane which 
is H| V{(j>) = 2f3 2 T 3 / cosh(V^m0) where (3 is 
a warp factor at the position of the D3 anti-D3 
in the internal compact space, T3 is the tension 
of branes and m is the string mass scale. At 
the late time (4> ~~ * 00 ) the potential behaves 
as V(<p) ~ /3 2 T 3 e~^ m ' i '. 

- (iiic) V(<P) = V e-^ M2 <t> 2 . 

This is the case in which the tachyon rolls 
down toward <j) — > 00 unlike the potential (iic) . 
A has a dependence A = M }) M 2 V r o ~ X/2 0e* Af2 * a 
thereby giving A — * 00 as <p ~ * 00 • 

It may be noted that the potentials (iiib) & (hie) 
listed above can be obtained in the frame work of 
string tachyons whereas the inverse power-law type 
potentials (iiia) are motivated from purely phe- 
nomenological considerations. 

We can now summarise the above behaviour in terms of 
conditions on V ((/)). The system approaches A — > when 
the slope of the potential is less steep than that of V = 
M 2 <f>~ 2 and the field 4> evolves toward infinity without 
oscillations [(iia) and (iib)]. The case (iic) corresponds 
to the one in which the field <f> oscillates as it evolves 
toward zero asymptotically. 

On the other hand, we have an asymptotic value 
|A| — > 00 when the potential is steeper than that of 
V = M 2 <\>~ 2 and the field (j) evolves toward infinity with- 
out oscillations [(iiia), (iiib) and (iiic)]. One may consider 
the potential V(<fi) = Vo((j> — 0o)™ j for n positive, that 
has a dependence A cx (0— 0o) - ' n ' 2+ j thereby showing 
a divergence of A for (j) — * 0q. However this potential 
has a number of problems such as a divergent negative 
mass m 2 s = (logV)^ as <fi — > 4>q, which leads to a vio- 
lent instability of perturbations in the context of tachyon 
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cosmology |35J . We do not regard this as a realistic dark 
energy tachyon potential. Note that if the potential has 
a positive constant energy at 4> = 0o then this system 
reduces to that of the case (iic). 



III. CONSTANT A 

Let us first consider the situation in which A is con- 
stant, i.e., case (i) in the previous section. The fixed 
points for this system can be obtained by setting x' — 
and y' = in Eqs. (jHJ and @. These are summarized 
in Table I. Essentially we have 4 fixed points: (a) x = 0, 
y = 0, (b) x = ±1, y = 0, (c) x = Xy s /V3, y = y s and 
(d) x = y/j, y = ±y/3j/\. Here y s is defined by 



VA 4 + 36 - A 2 



1/2 



(19) 



The cases (b) and (d) are divided into two cases, respec- 
tively, depending on the signs of x. The cases (c) and 
(d) correspond to stable fixed points. These satisfy the 

conditions 3a; = V3Xy and 37 = \/3Xxy - 1 - 3 ^ J ~ X ^ v 

Eqs. © and ©. 



A. Stability of the fixed point solutions 

We now study the stability around the critical points 
given in Table I. Consider small perturbations u and v 
about the points {x c ,y c ), i.e., 



x — x c + u , y = y c + v . 



(20) 



Substituting into Eqs. © and ©, leads to the first-order 
differential equations: 



= M 



(21) 



where M. is a matrix that depends upon x c and y c . 

The system can be regarded as being perturbatively 
stable when the the eigenvalues of the matrix M. are 
both negative :15|. In what follows we shall obtain the 
eigenvalues [i\ and fi2 for the fixed points in Table I and 
discuss their stability. 

• Case (a) [x c = 0, y c = 0]: 
The eigenvalues are 



/ii = -3 , [l 2 = 37/2 . 



(22) 



Therefore this critical point is an unstable saddle 
for 7 > 0, whereas it is a stable node for 7 = 0. This 
fixed point can not be used as a late-time attractor 
solution, since it leads to fi^, = 0. 



Case (b) [x c = ±1, y c = 0]: 
Since the eigenvalues are 

m = 6 , [i 2 = 37/2 . 



(23) 



this fixed point is an unstable node. This corre- 
sponds to a dust-like solution with 7^ = <p 2 = 1, 
but the system tends to repel from this critical 
point. 

• Case (c) [x c = Ay s /\/3, y c = y s \: 
The eigenvalues are 

A 2 



Mi 

/'2 



3+-(VA 4 + 36-A 2 ) 
A 2 



-37+ir(\/A 4 + 36-A 2 ) 



(24) 
(25) 



where [i\ ranges between — 3 < [i\ < —3/2. We 
have /12 < for 



7 > 7; 



A 2 



36 - A 



(26) 



This means that the fixed point is a stable node for 
7 > 7s, whereas it is an unstable saddle point for 
7 < 7s- 

• Case (d) [x c = ±^7, y c = ±y/3j/X\: 
The eigenvalues are 



Mi, 2 



7 - 2± W177 2 - 207 + 4- 



48 
A 2 



7 2 vT 



(27) 



The real parts of \i\ and p,i are both negative if the 
condition 



0< 7 <7, = ^(VA 4 + 36-A 2 ), 



(28) 



is satisfied. Note that 7 S is always smaller than 1. 
When the square root in Eq. H27JI is positive, the 
fixed point is a stable node. The fixed point is a 
stable spiral when the square root in Eq. I|27l) is 
negative. 

The values fl^ and 7^ at the critical point are 



r2<i, = 



30 



1 



A 2 % /T^ 



70 = 7, 



(29) 



which corresponds to a scaling solution in which 
the energy densities and pb decrease with the 
same rate. However we need to caution the reader 
in that the scaling solution does not exist in the 
cither the matter (7 = 1) or radiation dominated 
(7 = 4/3) eras, because the existence of the scaling 
solution requires the condition < 7 < 7s < 1- 
In this sense this solution can not be applied as a 
realistic model of dark energy. 
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19 


(a) 








All A and 7 


Unstable saddle for 7 > 
Stable node for 7 = 








fbl) 
\ UX J 


i 


o 


All A and 7 


Unstable node 


1 


1 


(b2) 


-l 





All A and 7 


Unstable node 


1 


1 


(c) 


\y s /V3 


y s 


All A and 7 


Stable node for 7 > 7 S 
Unstable saddle for 7 < 7s 


1 


A 2 J/ S 73 


(dl) 






A > and 7 < 7s 


Stable node 


37 1 

\' z \A-7 


7 


(d2) 


-V7 


-a/37/A 


A < and 7 < 7s 


Stable node 


37 1 

X ' 2 


7 



TABLE I: The critical points for constant A 



B. The existence of scaling solutions 

It was shown in Ref. |25| that the existence of scaling 
solutions restricts the form of the Lagrangian of a scalar 
field ip to be 



p(X,<p)=Xg(Xe X *), 



(30) 



where X = —g tll/ d IJ ,tpd u <p/2 and g is any function of 
Xe Xv . Equation JSDJl was derived by starting from a 
general Lagrangian p(X, ip) which is an arbitrary func- 
tion of X and ip. Recently this result was extended to a 
more general background given by H 2 cx p n [27| ■ 
The Lagrangian of our tachyon system is 



P(X, 



-V(<j))\/l-2X , 



(31) 



where X = — g^ d fJi (t>d v 4)/2. At first glance it seems that 
this lagrangian does not satisfy the condition for the ex- 
istence of scaling solutions given in Eq. 13Ufl . However we 
have earlier shown that this tachyon system l|31[) does in 
fact possess scaling solutions given by Eq. 

One can clarify the situation by expressing the La- 
grangian (|30|l in terms of the variable <f> = {2/ X) e Xip / 2 . 
Then Eq. ijSTlj) becomes 



(32) 



where f(X) = X g(X). The tachyon system l|31[) can be 
accommodated by choosing 



^)af 2 , f(X) cx Vl - IX. 



(33) 



Therefore scaling solutions exist for the Lagrangian l|31|) 
in the case of the inverse square potential, as of course we 
already knew. It is interesting to note that Eq. 1)30(1 com- 
pletely fixes the form of the Lagrangian for the existence 
of scaling solutions [2^, |2l} • 



In what follows, we would like to clarify several important 
points concerning the application of this model to dark 
energy. 

Employing slow-roll approximations: 3H(j) — —V^/V 
and 3MpH 2 ~ V(cf>) in a scalar-field dominated universe 
(p^ > pb), we obtain 



l 2M p 

a oc t p , with p 



M 



(34) 
(35) 



for the potential l|17|l . In order to have acceleration at 
late times, we clearly require p > 1, i.e., 



M > V2M V 



(36) 



Furthermore M needs to be much larger than the Planck 
mass to obtain a significant acceleration (p> 1). Such a 
large mass is problematic as we expect general relativity 
itself to break down in such a regime. This problem is 
fortunately alleviated for the inverse power-law potential 
V = M 4 ~ n cf)~ n with < n < 2, as we will see later. 

When n = 2 the parameter A defined in Eq. (|11() is a 
constant and given by 



M V P 



(37) 



Then we require A < \[2 for the acceleration. The fixed 
point (c) in Table I is the only stable attractor solution 
with an accelerating universe at late times (recall that 
the fixed point (d) is not a realistic solution). Let us 
consider the case of A <§; 1 (i.e., p S> 1). The fixed point 
(c) is approximately given as 



A 

x ~ —= , y ~ 
V3 

together with the equation of state 



12 



(38) 



Late time behavior 



The cosmological dynamics of the tachyon field with 
inverse square potential (|17|l was studied in Ref. 0, ^| . 



y_2_ 
70 ~ T ~ 3p~ 



(39) 



The condition for acceleration (p > 1) translates into 
70 < 2/3. One can easily verify that the slow-roll solution 
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in Eq. I|34(l is identical to the critical point (c) given in 
Eq. I|38l) . Therefore the slow-roll solution l|34[l is a stable 
attractor which gives Q$ — 1 and 7^ = 2/(3p). 

We note that there is another constraint coming from 
the energy scale of the tachyon potential. Since the en- 
ergy density of the tachyon is supposed to be the same 
order as the present critical density p cr — 10~ 47 GeV 4 , 
we require the condition M 2 ^ 2 ~ p cr . This restricts 
the present field values to be 4>qM > 10 60 . 



IV. CASE OF A 







In this section we shall study the case in which A dy- 
namically approaches 0. Examples of potentials which 
exhibit this behavior was presented in section II. In the 
limit of A — > Eqs. ijSJ and © read 



x' = — 3x(l — x 2 ) 

v =- lr y- 



(40) 
(41) 



By Eq. (gOl x> < for x > and x' > for x < (note 
that x 2 is in the range < x 2 < 1). Therefore x — > 0. 
Integrating Eq. I|41l) gives y cx l/H, which means that y 
continues to increase toward y — 1 as H decreases. Then 
the attractor solution should correspond to x — and 
y = 1. This discussion neglects the contribution from 
terms like ^/3Xy, but the solution actually approaches 
the attractor (x, y) = (0, 1) for A — > as we see below. 

While the fixed points we found in the previous section 
correspond to the case of constant A, we may regard these 
as "instantaneous" critical points as the function A(iV) 
evolves toward 0. The only stable critical point leading 
to an acceleration at late times is that of case (c) in Table 
I, where the dynamical critical point around A(iV) = is 
approximately given by Eq. I)38|l ■ This approaches x — > 
and y — ► 1 as A — > 0. Note also that 7^ — > as A — > 0, 
which corresponds to the equation of state of a cosmo- 
logical constant. We would like to numerically confirm 
that the system actually approaches the dynamical criti- 
cal point. In what follows we classify the situation in two 
classes: (A) A — > without any oscillations of <f> and (B) 
A — > with <j> oscillating. 



A 



with no oscillations of 



The potentials which lead to this behaviour correspond 
to (iia) V(<j>) = M 4 -"0-" with < n < 2 and (iib) 
V(<j>) = Voe 1/(M0) ■ In this class of models the field 4> ->■ 00 
without oscillations. As long as the tachyon potential is 
not steep relative to the inverse square potential, i.e., 
n < 2, X(N) asymptotically approaches as increases. 
This means that the equation of state of the field <fi ap- 
proaches 70 ~ X(N) 2 /3 — > 0, which leads to the universe 
accelerating at late times. 

Let us consider the inverse power-law potential V(<f>) — 
M 4 ~ n (j)~ n . In this case one has T = (n+l)/n in Eq. (fTT)) . 



which means that A continues to decrease for < n < 2. 
The slow-roll parameter for the tachyon-type scalar field 
is @ 



m p ( V A 1 _ n2 ( M v 



V V 



1 



M ) {(j)M) 2 



(42) 



When < n < 2, e decreases as the field evolves toward 
large values. The condition for the accelerated expansion 
corresponds to e < 1, which yields 



. „ , ( n M v 
V \V2M 



2/(2-n) 



(43) 



If this tachyon field is to be responsible for the observed 
inflation today, then it must satisfy Eq. ffify . 

The present potential energy is approximated as 
V(4> ) = M 4 /(<f> M) n ~ p CI ~ 10" 47 GeV 4 . Combining 
this relation with Eq. Q43|l we get 



M 



> 



Mi 



l-n/2 



V2 



1 1/(4-1) 



(44) 



The r.h.s. becomes smaller as n decreases. For exam- 
ple one has M/M p > 1(T 20 for n = 1. Therefore the 
super-Planckian problem for the inverse square potential 
is alleviated for n < 2. 

In Fig. n we plot the cosmological evolution for n = 1 
with initial conditions Xi = 0.8, g/j = 5.0 x 10~ 4 and 
Aj = 1.0. We consider a pressureless dust (7 = 1) as a 
background fluid. Note that e and A satisfy the relation 
e = A 2 /2 by Eqs. (Q3J and Our choice h = I 

corresponds to the initial condition = 0.5. This does 
not mean that inflation occurs at the initial stage, since 
the energy density of the barotropic fluid dominates over 
that of the scalar field. When the energy density of <f> 
eventually wins out, this leads to the acceleration of the 
universe since the slow-roll parameter is smaller than of 
order 1. 

From the panel (A) of Fig.^ we find that x decreases 
initially. This comes from the fact that the condition 
3x 3> V3\y holds in Eq. (JHJ for our initial conditions. 
On the other hand y grows through the relation y' « 
(3j/2)y. Since this growth is rather rapid and A is nearly 
constant initially, the v^Ay term temporarily surpasses 
the 3a; term in Eq. <(SJ), which leads to the increase of x 
for a short period. After that, A begins to decrease [see 
the panel (C)] and this has the effect of balancing the two 
terms (3x » ^/3\y). In Fig.Q]we plot the dynamically 
changing critical points (c) in Table I, i.e., (x c ,y c ) — 
(X(N)y s (N)/V3,y s (N)) with y s (N) = [(^A(7V) 4 + 36- 
A(iV) 2 )/6] 1 / 2 . We find that the solution approaches these 
"instantaneous" critical points with the decrease of x. 
Therefore the discussion of constant A can be applied to 
the case of varying A after the system approaches the 
stable attractor solutions. One can find the asymptotic 
evolution of A by substituting Eq. for Eq. (fTUfl . This 
gives the dependences A cx 1/%/a", x cx 1/s/N and 1 — y cx 
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1/N 2 , which we confirmed numerically. Note that this 
relation holds as long as T is asymptotically constant 
with r < 3/2. 

The evolution of fi^, Qg, A, 7^ and q is plotted as well 
in the panels (B) and (C) in Fig. ^ The deceleration 
parameter q becomes negative for N > 100, after which 
the system enters the acceleration stage with the growth 
of f20. We checked that 74, evolves toward with the 
decrease of X(N) keeping the relation 7^ ~ X(N) 2 /3. 

We can also account for a combined system of two flu- 
ids (matter and radiation) with a scalar field <j>. We have 
done this, running our numerical code from a redshift 
2 = 10 6 , finding that it is possible to obtain a viable cos- 
mological evolution for the inverse power-law potential 
with n < 2. The basic property of the dynamical system 
for the potential V(<p) = V^e 1 ^^ is similar to what we 
discussed above. 
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B. A — > with oscillations of (f> 

One example exhibiting this type of behaviour is the 
case (iic), i.e., a rolling massive scalar field with potential 
V((f) = Vqc^ m $ |2^|. This potential has a minimum at 
4> = with an energy density Vq- As an example, consider 
an anti-D3 brane at the tip of Klebanov-Strassler (KS) 
throat in the KKLT setup Q • Since there is a warp factor 
f3 at this point, the potential of a massive excitation of the 
anti-D3-brane (f> may be written as V(</>) = [3 2 T^,e^ l3m $ , 
where T3 is the brane tension and m is the mass of the 
excited state of the brane which is of order the string 
mass scale. Hence, if we introduce a very small warp 
factor (3, it should be possible to explain the origin of 
the present dark energy. Note that with a warp factor f3 
of order 1, the massive scalar field decays to <\> = very 
soon in the reheating epoch. However, for a very small 
warp factor, the stabilization of the field <j> is considerably 
delayed allowing it to play the role of dark energy today. 

Another possibility is to consider an anti-D3-brane 
with a warp factor (3 1 and a negative cosmological con- 
stant (—A) arising from the stabilization of modulus 
fields in the KKLT vacua with the assumption that the 
potential energy of the rolling massive scalar does not ex- 
actly cancel the cosmological constant (V(4>) — A ~ p cx ) 
|29j . The rolling massive scalar, in order to play the roll 
of quintessence, is expected to be stabilized in future. 

In this work we shall simply adopt the potential 
V{4>) = Voe2 M v and study the dynamics of the sys- 
tem when the field oscillates around the potential mini- 
mum. Since A = — M p M 2 0/v / Voe _ 3 M * , this quantity 
gradually decreases toward and is accompanied with 
the oscillations of cf>. In Fig.|3 we plot one example of 
the cosmological evolution for this scalar potential. The 
quantity x approaches with damped oscillations as ex- 
pected. We find that the "instantaneous" critical points 
(c) in Table I provide a good description of the late time 
evolution of x and y. fl^ begins to grow toward 1 around 




-0.20 I J 

1 100 1 J 1 6 10 s 1 1 

N 

(C) 

1 .2 I i 




-1.2 I J 

1 100 1 J 1 ft 1 8 1 1 
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FIG. 1: Cosmological dynamics for the field <j> with potential 
V((f>) = M 3 cf>~ 1 and a barotropic fluid with 7=1. We choose 
initial conditions as Xi = 0.8, yt = 5.0 x 10~ 4 and Xi = 1.0. 
Each panel corresponds to the evolution with (A) x, y and 
the critical points x c , y c for the case (c) in Table I, (B) Q^, 
JIb and (C) A, 7^ and the deceleration parameter q. 
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N ~ 20, after which the system enters the accelerating 
stage [see the evolution of q in the panel (C)]. 



V. CASE OF |A| oo 

There are a number of potentials that give | A| — ^ oo 
asymptotically. One can classify the tachyon potentials 
into two classes: (1) |A| — > oo without the oscillation of 
the field <j> [(iiia), (iiib), (hie)] and (2) |A| — > oo with the 
oscillation of the field, that is, V(<j)) — Vq((/) — <fio) n for 
positive n. In the latter case the effective mass squared 
ml s = (log V)^ of the field <f) is negative, which leads to 
a violent instability for the tachyon perturbations |35| . 
This is not regarded as a stable attractor unlike the po- 
tential V(4>) — Vbe2 M v discussed in the previous sec- 
tion. 

We shall investigate the case (1) in which |A| continues 
to grow asymptotically. For example let us consider the 
exponential potential V(4>) = Vbe _M * with > [case 
(iiib)]. Since T — 1 for this potential, we have A' = 
(\ // 3/2)X 2 xy, thereby leading to the growth of A for x > 
0. In Fig.[3]we plot the cosmological evolution for this 
system with initial conditions Xi — 0.7, yi = 1.0 x 10~ 4 
and Xi — 0.5. There is a short initial stage in which the 
conditions, x <C 1 and y <C 1, are satisfied with ~ 
and 70 ~ 0. This is an unstable fixed point (a) in Table I, 
thus showing a deviation from x = and y = for N > 
100. From the panel (C) we find that the universe begins 
to accelerate for N > 400 (see the evolution of q) . After 
that x and y approach the instantaneous critical points 
(c) in Table I. As long as A < 0(1), one can see that the 
condition q < is satisfied by substituting x s = Xy s /\/2> 
and y = y s for Eq. l(TB|) . With the growth of A, however, 
the acceleration of the universe stops for N > 6.5 x 10 4 . 
At this stage x grows toward 1, where as y decreases 
toward 0. The asymptotic behavior corresponds to x — > 1 
and y — > V3/A, which can be obtained by taking the 
limit A ^> 1 for the instantaneous critical points (c) in 
Table I. Substituting this solution for Eq. IjlOII . we obtain 
A oc e 3Ar / 2 . This relation holds when T is asymptotically 
constant with T > 3/2. In this regime the equation of 
state for the field 4> is characterized by 7^ ~ 1 (i.e., a 
pressureless dust) with a dominant energy density (f^ ~ 
1)- 

It is worth reflecting here briefly on the fact that there 
exists a period of transient acceleration when we have 
an exponential potential. Because of a dynamical change 
of A, it is possible to have a temporal acceleration for 
A < 1 and have a deceleration for A 3> 1. If this tem- 
poral acceleration corresponds to the one at present, the 
universe will eventually enter the non-accelerating regime 
in which the tachyon field behaves as a pressureless dust. 
This is a nice feature of the model which makes it free 
from the future event horizon problem present in most of 
the quintessence models. 

As long as |A| approaches infinity asymptotically, we 
do not have an acceleration of the universe at late times. 
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FIG. 2: Cosmological dynamics for the field <j> with potential 
V{4>) = Voe3 M2|/>2 and a barotropic fluid with 7=1. The 
initial conditions are Xi = —0.5, yt — 6.7 x 1CF 3 and A; = 
—7.8. We show the evolution of the same quantities as in 
Fig. 1. 
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FIG. 3: Cosmological evolution for the field (j> with potential 
V(4>) — Vbe _A "^ and the barotropic fluid with 7 = 1. Initial 
conditions are Xi = 0.7, yi = 1.0 x 10~ 4 and Ai = 0.5. We 
show the evolution of the same quantities as in Fig. 1. 



This comes from the fact that the fixed point (c) in Ta- 
ble I approaches x — > 1 and y ~ \/3/A — > as |A| — > 00. 
Therefore any potentials which exhibit this behavior can 
not be used as a late-time dark energy candidate. In 
order to lead to acceleration at late times, we require 
that the tachyon potential has the property (i) A = const 
or (ii) A — > asymptotically, as we showed in previous 
sections. Of course, as we have just mentioned, it is pos- 
sible that the acceleration we are experiencing is not a 
late time evolution, but part of a transient regime! 



VI. SUMMARY 

In this paper we have studied the cosmology associ- 
ated with the tachyon field (f> as the dark energy. In the 
tachyon system the inverse square potential l|17fl is the 
marginal case which leads to an accelerating universe. 
However, it is not the only possible tachyon potential 
and so we have deliberately set out to develop a unified 
analysis that can be applied for any type of tachyon po- 
tentials. 

A crucial quantity to determine the cosmological dy- 
namics is A defined in Eq. (|llfl . The inverse square poten- 
tial l|17fl corresponds to a constant A with A = 2M p /M. 
In this case we have essentially four critical points as 
shown in Table I. Among them a viable fixed point that 
gives a stable attractor for dark energy is the case (c) 
with a barotropic equation of state 7 larger than "f s [j s 
is defined in Eq. 1)26(1] . The critical point (d) corresponds 
to a scaling solution in which the energy density of <fi de- 
creases similarly to that of the barotropic fluid. However 
this is not a viable cosmological scaling solution, since 
the existence of it requires the condition 7 < 1 . 

The tachyon potentials can be placed into three classes: 
(i) A = const, (ii) A — > asymptotically and (iii) |A| — > 00 
asymptotically. The class (ii) corresponds to the case 
in which the potential is not steep relative to V(<f>) — 
AI 2 (f)~ 2 , whereas the potential in the class (iii) is steeper 
than V{4>) = M 2 cj)~ 2 . While the former leads to the 
acceleration at late times, the latter does not. 

We have carried out a detailed analysis about the cos- 
mological evolution of the tachyon system in two cases: 
(1) A — > without the oscillation of <f> and (2) A — * 
with the oscillation of <f>. We adopt the inverse power- 
law potential V{4>) = M A ~ n 4>~ n with < n < 2 as an 
example of the case (1). This is favorable relative to 
the inverse quadratic potential [n = 2), since the mass 
scale M is not severely constrained. We solved the dy- 
namical equations (|5|)- (|l(Jfl numerically and found that 
the solutions approach the "instantaneous" critical point 
(c) in Table I with the decrease of A toward [see the 
panel (A) of Fig. 1]. Since x ~ A(A)/ % /3 -> and 
y ~ 1 - A(A) 2 /12 -t las X(N) -> 0, the universe 
exhibits an acceleration at late times. This is a dark 
energy scenario in which the future universe is domi- 
nated by the field <f> (f2^ ~ 1) with an equation of state 
70 ~ 0. We adopt the rolling massive scalar potential 
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V{(/>) = V Q e2 M <t> to study the case (2) mentioned above. 
We find that x and y again approach the instantaneous 
critical point (c) with oscillations. Since the potential 
has an energy Vo at the potential minimum, this eventu- 
ally leads to an acceleration of the universe even if the 
scalar field oscillates. 

On the other hand we do not find an accelerated ex- 
pansion at late times when |A| grows towards infinity, 
since the dynamical critical point (c) approaches x — > 1 
and y — > as |A| — > oo. Nevertheless we found that a 
transient acceleration can occur in the region where X(N) 
is smaller than of order unity (see Fig. 3). It is tempt- 
ing to speculate that this could be the transient regime 
we are experiencing today, and to try and predict when 
we expect it to change over again to a matter dominated 
era. Recall that the supernova data is not actually in- 
forming us about the evolution today (z = 0), rather it 
tells us that the universe is accelerating at a redshift of 
order 0.05. Perhaps we are actually back in a matter 
dominated regime today and just don't know it yet! 

Inspite of the existing features of cosmological dynam- 
ics based upon DBI scalars, it may happen that these 
models lead to the formation of caustics where the second 
and the higher-order derivatives of the field become sin- 
gular. We do not know whether caustics are generic pre- 
diction of string theory or appear as a result of the deriva- 
tive truncation leading to the DBI action. As demon- 
strated in Ref. |36|. caustics inevitably form in tachyon 



system with potentials decaying as <fi~ 2 or faster at in- 
finity. It remains to extend the analysis of Ref. j3|J to 
the case of potentials such as V ~ <fi~ n with < n < 2 
analysed in our paper. Caustics normally form in systems 
with pressureless dust which is mimicked by tachyon field 
with run away potentials. It is therefore quite likely that 
caustics may not develop in Born-Infeld systems with a 
ground state at at a finite value of the field. The rolling 
massive scalar potential V((f>) — Vbe3 M ^ belongs to 
this category. In our opinion, these are important issues 
which require further investigations. 

Finally we believe it is quite intriguing that the tachyon 
system provides rich and fruitful cosmological scenarios 
for dark energy. The classification we have performed in 
this paper provides a very useful way to find out about 
the cosmological evolution for any type of tachyon po- 
tentials. 
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